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I.  INTRODUCTION 

In  the  last  couple  of  years  the  search  for  new  materials  that  possess 
unique  physical  properties  has  been  going  through  a  revival  phase.  Especially 
the  development  of  novel  sophisticated  experimental  techniques,  such  as  laser 
desorption  (LO)  and  chemical-vapor  deposition  (CVO)  have  led  to  the 
generation  of  hitherto  unknown,  unusual  compounds.  Amongst  these  materials, 
carbon-  and  silicon-containing  nitrides  have  received  special  interest.^  Thus, 
the  covalently  bound  polymer  carbon(iV)  nitride  has  only  recently  been 
synthesized  by  a  combination  of  LD  and  CVD  techniques.^  This  structurally  not 
yet  fully  characterized  substance  seems  to  exhibit  the  theoretically  predicted 
physical  properties  such  as  an  extreme  hardness  and  diatomic-like  modules  of 
elasticity.  It  is  believed  to  possess  a  structure  similar  to  the  binary  silicon  nitride 
p-Si3N4. 

In  general,  it  is  known  that  the  p-block  main-group  elements  form  nitrides 
with  exceptional  chemical  inertness  and  mechanical  hardness.^  However, 
binary  silicon  nitrides  are  not  only  of  interest  to  material  scientists.  Small  silicon- 
containing  clusters  are  believed  to  be  formed  in  silicon-etching  processes  and 
thus,  as  they  can  also  be  viewed  upon  as  prototypes  for  the  interaction  of  main- 
group  elements  with  silicon  sites  in  clusters  and  surfaces,  have  received  much 
attention  by  experimentalists^  as  well  as  theoreticians.^ 

To  our  knowledge  there  exists  only  very  little  information  on  small  binary 
silicon-nitrogen  clusters.  Recently,  Brodie  et  al.^  have  carried  out  pyrolysis- 
solid-state  NMR  studies  on  silazanes  in  order  to  obtain  mechanistic  insight  into 
the  genesis  of  the  siliconnitride  solid  (P-Si3N4),  and  some  of  us  have  been  able 
to  characterize  the  Si2N  and  Si20  cluster  molecules^  by  neutralization- 
reionization  mass  spectrometry  (NRMS).®  This  mass  spectrometric  technique 
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has  been  applied  successfully  for  the  generation  of  numerous  elusive  silicon- 
containing  molecules,  which  due  to  the  facile  intermolecular  reactions  could  not 
be  studied  in  the  condensed  phased 

In  this  publication  we  wish  to  report  the  results  of  combined  experimental 
and  ab  initio  investigations  on  the  hitherto  unknown  SisN  cluster  and  its  cationic 
oi3N+  and  anionic  SiaN-  counterparts.  The  cationic  and  neutral  species  have 
been  structurally  characterized  by  collisional-activation  and  neutralization- 
reionization  mass  spectrometry.^  High-level  ab  initio  methods  have  been 
employed  for  the  structural  and  energetic  description  of  the  various 
experimentally  observed  neutral,  cationic  and  anionic  SipN  and  Sin  (n=1-3) 
molecules.  Calculated  structural  data  of  the  experimentally  observed  species 
will  be  given  together  with  vibrational  spectra  to  facilitate  the  spectroscopic 
observation  of  these  novel  molecules. 

II.  EXPERIMENTAL  DETAILS 

The  collision  experiments  were  performed  in  a  substantially  modified 
ZAB  four  sector  tandem  mass  spectrometer  with  BEBE  modification  (B  stands 
for  magnetic  and  E  for  electrostatic  analyzer).  For  a  detailed  description  of  the 
mass  spectrometer  the  reader  is  refereed  to  previous  publications.^  For  the 
cluster  generation,  trisilylamine,  N(SiH3)3  was  used  as  the  precursor  molecule. 
This  compound  was  prepared  via  standard  literature  procedures.  Care  had  to 
be  taken  to  avoid  possible  leaks  in  the  inlet  system  because  oxygen  traces  can 
lead  to  sudden  explosion  of  the  amine.  The  sample  was  introduced  into  the  ion 
source  of  the  mass  spectrometer  through  a  chemical  ionization  gas  inlet  system. 
The  cationic  cluster  could  easily  be  obtained  upon  dissociative  electron  impact 
ionization  of  N(SiH3)3.  The  source  conditions  were  as  follows:  electron  energy 
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50  eV;  repeller  voltage  ca.  30  V;  ion  source  temperature,  2000C;  accelerating 
voltage  8  kV.  The  ions  of  interest  were  mass  selected  by  means  of  (B1)(E1). 
Collision  experiments  were  performed  with  80%  transmittance  (T)  of  the  ion 
beam  in  order  to  ensure  single-collision  conditions.  In  the  case  of  the  collisional 
activation  (CA)  mass  spectrum,  helium  was  introduced  as  a  collision  gas  into 
the  second  part  of  a  differentially  pumped  tandem  collision  cell.  For  the 
neutralization-reionization  experiments  xenon  was  used  as  the  neutralization 
agent  in  the  first  part  of  the  tandem  collision  cell.  The  remaining  ions  were 
electrostatically  deflected  from  the  beam  by  means  of  an  electrode  (charged  to 
1 000  V)  situated  between  the  two  collision  cells.  The  neutral  molecules  were 
reionized  by  colliding  the  beam  with  oxygen,  For  a  particle  m/z=98  having  been 
accelerated  to  8  keV-energy  it  takes  approximately  6  microseconds  to  travel 
from  the  first  to  the  second  collision  ceil,  thus  providing  a  lower  limit  for  the 
minimal  lifetime  of  the  neutral  molecule.  The  ionic  dissociation  products  were 
recorded  by  scanning  the  second  magnetic  analyzer  (B2).  To  improve  the 
signal-to-noise-ratio  of  the  spectra  the  data  of  20-50  accumulated  scans  were 
on-line  processed  with  an  AMO-Intectra  data  system. 

11.  COMPUTATIONAL  DETAILS 

The  bond  lengths  and  angles  of  neutral,  cationic  and  anionic  SiN,  SiaN, 
and  SiaN,  as  well  as  Sia  and  Sia  were  optimized  by  employing  analytical 
gradients^  ^  using  the  Gaussian  92  program^^  with  polarized  split-valence  basis 
sets  (6-31 at  the  SCF  and  the  MP2(full)  level  (USCF  and  UMP2(full) 
for  open  shell  systems  of  theory).  Our  optimal-geometric  parameters  are 
presented  in  Fig.  1  and  our  simulated  IR-spectra  based  on  the  MP2(full)/6- 
31 1 +G*  calculated  frequencies  and  their  intensities  are  drawn  in  Fig.  2. 
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The  MP2(full)/6-31 1 +G*  equilibrium  geometries  were  used  to  evaluate 
electron  correlation  corrections  in  the  frozen-core  approximation  by  full  fourth 
order^s  Moller-Plesset  perturbation  theory  and  by  the  (U)QCISD(T)  method^® 
using  6-311+G(2df)  basis  sets.  The  USCF  wave  functions  for  open  shell 
systems  were  projected  to  pure  spectroscopic  states  for  which  the 
corresponding  results  are  denoted  PUSCF,  PMP2,  PMP3  and  PMP4.17 
Calculated  reaction  energies  and  available  experimental  data  are  presented  in 
Table  I. 

IV.  THEORETICAL  RESULTS 

SIN,  SIN-**  and  SINvThe  neutral  SiN  molecule  is  known  to  have  a 
doublet  2£+  (lo22c2l7c^3oi)  ground  electronic  state.^8-19  The  lowest  excited 
state  2ni  (1o22o21;c33o2)  is  only  0.25  eV  higher  in  energy.''^  At  our  highest 
level  (QCISD(T)+ZPE),  the  2l+  state  is  indeed  found  to  be  the  ground  state  and 
the  first  excited  2ni  state  is  0.35  eV  higher  in  energy.  Our  calculated  bond  length 
Re(Si-N)=1 .53  A  for  the  ground  2^+  state  at  the  MP2(full)/6-31 1 +G*  level  is 
somewhat  shorter  than  the  experimental  value  1.57  A.i®-^^  To  the  best  of  our 
knowledge,  experimental  structural  data  on  SiN+  and  SiN-  are  not  known. 

From  the  2^+  and  2ni  electronic  states  of  the  SiN  molecule,  the  two 
cationic  triplet  (1o22a2l7c23o2)  and  ^II  (1a22a2l7t33ai)  states  are  expected 
to  be  the  most  stable.  Bruna,  Peyerimhoff  and  Buenker  (BPB)20  studied  five 
low-lying 

(1o22o2iji4)  and  U  (1o22o2l7c23o2)  electronic  states  of  SiN+  using  large-scale 
Cl  calculations  and  found  that  the  is  the  ntost  stable  and  the  first  state  is 
0.21  eV  higher  in  energy.20  We  repeated  the  calculations  of  the  both  the 
and  3n  states  at  the  MP2(full)/6-31 1+G*  level  and  our  bond  lengths  Re(Si- 
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N)»1.78  A  (3I-)  and  R0(Si-N)=1 .65  A  (^n)  reasonably  agree  with  the  Re(Si- 
N)=1.81  A  (3S-)  and  Re(Si-N)=1.71  A  (3n)  BPB  data.^3  Our  (3l+  3n) 

excitation  energy  (0.13  eV)  is  also  close  to  the  0.21  eV  BPB  value.  Our  first 
adiabatic  IE  for  SiN  is  10.0  eV  at  the  QCISD(T)  level. 

The  experimental  electron  affinity  is  not  known  for  SiN.  The  lowest 
singlet  anion  electronic  state  (lo^aa^ijc^Sa^)  may  be  derived  from  either 
the  ground  or  the  first  excited  electronic  state  of  neutral  SiN  when  an  electron  is 
added  to  the  lowest  energy  available  MO.  The  isoelectronic  SiO  molecule  has  a 
singlet  ground  state  with  its  first  excited  state  lying  more 

than  4  eV  higher  in  energy.iQ  Therefore  we  studied  only  the  ^1+  state  for  SiN'. 

The  optimized  bond  lengths  and  vibrational  frequencies  at  the 
MP2(full)y6-31 1 +G*  level  are:  1.63  A  and  1005  cm*^  which  agree  well  with 
results  of  Peterson  and  Woods^i  Re(Si-N)= 1.605  A  and  o)e*1092  cm'i  for  SiN* 
at  the  MP4  level.  The  first  adiabatic  electron  affinity  of  SiN  is  found  to  be 
3.36  eV  at  the  PMP4  level  and  2.91  eV  at  the  QCISD{T)  level,  both  of  which  are 
lower  than  the  experimental  EA  (3.821  eV)22  of  the  valence  isoelectronic  CN 
molecule.  The  electron  affinity  of  SiN  calculated  indirectly  by  Peterson  and 
Woods  is  3.32  eV.2i  Dissociation  energies  (at  the  QCISD(T)level)  for  SiN-^  into 
Si+  +  N,  SiN  into  Si  +  N  and  SiN*  into  Si*  +  N  are:  52  kcal/mol,  96  kcal/mol  and 
1 35  kcal/mol,  respectively. 

Si2N,  Si2N'*’  and  Sl2N'.  The  neutral  Si2N  molecule  and  SiaN-*-  cation 
were  recently  observed  in  mass  spectrometric  studies,^®  but  to  our  knowledge, 
there  exist  no  analogous  experimental  data  on  SiaN*,  and  theoretical  data  are 
not  available  for  any  of  these  species. 

We  optimized  the  geometry  for  the  cation  Si2N+  using  the  MP2(full)/6- 
31 1+G*  method  starting  from  an  angular  Cav  (^Ai)  structure  with  an  initial  SiNSi 
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angle  equal  to  140°  We  converged  to  Re(Si-N)=l  .631  A,  and  a  valence  angle 
of  Si-N-Si=180o.  That  is,  the  optimal  structure  I  (  geometry  shown  in  Fig.  1  and 
simulated  IR-spectrum  in  Fig.  2)  of  SiNSi+  is  linear  Dooh  : 
1ag2lou^2og2ljtu^2ou2l7Cg0),  which  is  different  from  the  isoelectronic  SiCSi 
molecule,23.24  which  has  a  singlet  bent  C2v  (^Ai)  structure.  Our  most  favorable 
dissociation  energies  (at  the  QCISD(T)  level,  see  Table  1)  for  Si2N+  is  into  Si-*-  + 
SiN  (+167  kcal/mol). 

We  optimized  the  geometry  for  both  linear  Dooh  (^Flg:  log^lou^ 
2ag2lnu^2au2l7Cgi)  and  bent  C2v  (^Ai:  1ai2lb222ai2ibi2 3ai22bi24aii) 
structures  of  Si2N.  At  the  MP2(full)/6-31 1+G*  level,  the  C2v  is  a  minimum  (see 
the  structure  II  in  Fig.  1  and  the  simulated  IR-spectrum  in  Fig.  2),  but  the  Dooh 
(2ng)  structure  III  has  one  imaginary  frequency  whose  eigenvector  leads  to  a 
linear  Coov  (^O  :  structure  IV  which  is  only  0.36  kcal/mol 

lower  in  energy  at  this  level  of  theory.  Optimization  of  the  Coov  (^II)  structure 
within  the  MP2(full)/6-31 1 +G*  framework  gives  a  minimum  with  Iwa 
nonequivalent  bond  lengths  RofSii -Nisi  .573  A  and  Re(Si2-N)=1 .694  A.  The 
distortion  of  S2N  from  Dooh  into  Coov  symmetry  at  MP2(full)  level  is  due  to  an 
artifactual  symmetry  breaking.  For  our  large  6-31 1+G(2df)  basis  sets,  the  energy 
of  the  linear  Dooh  structure  III  is  lower  at  all  correlated  levels  (from  PMP2  to 
QCISD(T))  than  the  C«,v  structure  IV  using  MP2(full)/6-31 1+G*  geometries  for 
both  stnjctures.  Therefore  we  optimized  the  geometry  of  linear  SiNSi  at  the 
QCISD(T)/6-31 1 +G*  level,  where  we  found  no  symmetry  breaking  problem.  The 
linear  Dooh  (^Ilg)  structure  is  a  minimum  with  Re(Si-N)=1.643  A  and  vi(Og)=613 
cm-1,  v2(jcu)=221  cm-^  v2'(Jtu)»130  cm*^  ano  V3(au)=898  cm-i  (simulated  IR- 
spectrum  is  shown  in  Fig.  2).  Therefore  the  global  minimum  of  the  Si2N 
molecule  has  a  linear  Dooh  (^Ilg)  structure. 
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Both  II  and  III  structures  of  SiaN  are  very  close  in  energy  and  may  coexist 
in  the  gas  phase.  According  to  our  calculations,  the  energetically  favored 
dissociation  path  for  SiaN  is  into  SiN  +  Si  (+125  kcal/mol);  dissociation  into  Sia 
+  N  (+152  kcal/mol)  is  less  favorable.  Finally,  the  adiabatic  ionization  energy  of 
SiaN  is  predicted  to  be  6.2  eV  (see  Table  1). 

The  SiaN'  anion  is  isoelectronic  with  the  neutral  SiaO  molecule  and, 
according  to  previous  calculations,^^  the  angular  Cav  (^Ai)  structure  of  SiaO  is 
the  most  stable,  but  two  linear  SiOSi  (Dooh  and  SiSiO  (Coov  ,^I*) 
structures  are  18-27  kcal/mol  higher  in  energy.  We  optimized  the  geometry  for 
ail  of  these  three  structures  of  SiaN':  V  (Cav .  ^Ai:  1ai2iba22ai2lbi23ai2 
2bi24ai2).  VI  (Dch.  3lg-:  1ag2lOu22og2lKu42ou2l7Cg2)  and  VII  (Coov.  : 
1a22o23a2l7c^4a22jc2),  and  found  that  all  three  structures  are  local  minima  with 
the  singlet  V  structure  being  the  most  stable,  as  in  the  case  of  SiaO.^^  However, 
the  linear  VI  structure  is  only  5.0  kcal/mol  higher  In  energy  at  the  PMP4/6- 
31 1  +G(2df)  level  and  1 .7  kcal/mol  at  the  QCISD(T)/6-31 1  +G(2df)  level.  The  VII 
structure  is  39.6  kcal/mol  (at  the  QCISD(T)/6-311+G(2df)  level)  higher  in  energy 
than  the  most  stable  V.  The  geometries  of  all  these  structures  are  presented  in 
Fig.  1  and  the  simulated  IR-spectra  are  shown  in  Fig.  2. 

These  energy  data  indicate  that  the  SiaN*  anion  may  coexist  in  two 
molecular  forms  in  the  gas  phase.  Electron  photodetachment  experiments  of 
SiaN*  may  allow  experimental  verification  of  these  predictions.  Our  calculated 
adiabatic  electron  affinity  is  EAa(SiaN)=1.1  eV  for  SiaN.  Our  most  favorable 
dissociation  energies  (at  the  QCISD(T)  level,  see  Table  1)  for  SiaN+  into  Si+  + 
SiN,  SiaN  into  Si  +  SiN  and  SiaN*  into  SIN*  +  Si  are:  167  kcal/mol,  125  kcal/mol 
and  80  kcal/mol,  respectively. 
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SiaN,  SiaN-^  and  SiaN‘.  We  initially  performed  calculations  on  SiaN-*- 
at  the  SCF/6-31 UG*  level  and  found  the  high  symmetry  Dan  structure  (with 
Iai‘2le'^2ai’2la2‘'22e"*3e'2  electronic  configuration)  to  not  be  a  local 
minimum  either  in  its  triplet,  or  singlet  states.  The  triplet  Dan  (^A')  structure  VIII 
(Fig.  1;  nas  two  imaginary  frequencies  with  two  eigenvectors  leading  to  the  Y- 
type  structure  IX  and  to  the  T-type  structure  X  (both  have  the  same  electronic 
configuration:  Cav.  ^82,  (1ai2lb2^2ai23ai2lbi24ai22b2^  3b2^5ai^).  The 
structures  IX  and  X  (Fig.  1 )  are  local  minimum  and  saddle  point,  respectively. 

The  barrier  for  intramolecular  preudorotation  (top  of  the  barrier  corresponding  to 
the  structure  X)  is  3.4  kcal/mol  at  the  SCF/6-31 1  +G*+ZPE  level.  The  other  triplet 
(C2v.  ^Bi,  Iai2ib222ai23ai2ibi24ai22b22  2bi’5aii)  state  is  a  local  minimum 
at  the  Y-type  structure  XI.  This  structure,  relative  to  the  IX  structure,  is  2.4 
kcal/mol  (at  SCF/6-31 1 +G*)  higher  in  energy.  However  both  of  the  triplet 
minima  are  higher  in  energy  than  the  singlet  surface  of  the  SiaN-*-  cation. 

The  high  symmetry  singlet  Dan  (^Ai  with  1ai’2le'^2ai'2la2"22e’^3e'2 
electronic  configuration)  structure  is  not  expected  to  be  a  minimum  because  of 
Jahn-Teller  distortion  which  leads  to  the  Y-type  (XII)  and  the  T-type  (XIII)  of 
structures.  The  singlet  XII  structure  (Cav.  1ai22ai2lb223ai2ibi24ai22b22 
5ai2)  is  found  to  be  a  saddle  point,  whereas  the  singlet  XIII  structure  (Cav. 
1a^2ib222ai2ibi23ai22b224ai23b22)  is  a  global  minimum  for  SiaN+  at  the 
SCF/6-31 1 +G*  level.  The  barrier  for  intramolecular  pseudorotation  is  7.0 
kcal/mol  (SCF/6-31 1+G*+ZPE).  At  the  SCF  level,  the  singlet  potential  energy 
surface  is  lower  in  energy  than  the  lowest  triplet  302  potential  energy  surface  by 
more  than  10  kcal/mol.  Because  the  singlet  surface  is  expected  to  be  even  more 
stable  at  correlated  levels,  we  studied  only  the- singlet  structures  XII  and  XIII  at 
the  MP2(full)/6-311+G*  level. 


We  optimized  the  geometry  and  calculated  the  harmonic  frequencies  for 
the  singlet  XII  and  XIII  structures  at  the  MP2(full)/6*31 1 +G*  level.  These 
structures  are  shown  in  Fig.  1  and  the  vibrational  spectrum  of  XIII  is  shown  in 
Fig.  2.  At  MP2(full)/6-31 1+G*  level,  XIII  is  a  global  minimum  (similar  to  the 
isoelectronic  SiaC  molecule25.26)  and  XII  is  a  saddle  point.  However  at  this 
level  the  pseudorotation  barrier  is  substantial  (17.3  kcal/mol  including  ZPE 
correction).  The  global  minimum  structure  XIII  was  calculated  at  the  QCISD(T) 
and  PMP4  levels  to  obtain  a  more  accurate  dissociation  energy.  According  to 
our  calculations,  among  all  possible  dissociational  channels  SiaN-*-  into  Si  -«■ 
Si2N+  (66  kcal/mol  at  QCISD(T)  see  Table  1)  has  the  lowest  energy. 

From  the  order  of  the  orbitals  in  SiaN+  (Dah),  we  predict  that  the  neutral 
SiaN  molecule  has  occupation  1ai’2le"^2ai*2la2''22e'^3e'3  in  the  high 
symmetry  Dan  (^E' )  structure,  which  is  not  expected  to  be  a  minimum  because 
of  the  Jahn-Teller  distortion,  which  leads  to  the  Y-type  (202, 
Iai22ai2ib223ai2ibi24ai22b225ai2  3b2^2bi0la2°)  and  T-type  (2Ai, 

1  ai  22ai  2 1  b223ai  2 1  bi  24ai 22b225ai  ^ 3b22  2bi  o  1 32°)  structures.  Because 
orbitals  other  than  the  5ai  and  3b2  low-lying  orbital  may  play  important  roles, 
we  further  studied  ail  configurations  obtained  by  distributing  three  electrons 
among  the  5ai,  3b2,  2bi,  and  1a2  low-lying  orbitals  with  doublet  or  quartet 
coupling  of  spins. 

The  doublet  coupling  leads  to  twelve  states:  C2v  symmetry  2Ai 
(5aii3b222bi°1a2°,  structure  XIV),  2Ai  (5ai^3b2°2bi2ia2°,  structure  XV),  2Ai 
(5aii3b2°2bi°1a22,  structure  XVI),  202  (5ai23b2^2bi°1a2°,  structure  XVII),  202 
(5ai03b2^2bi2ia2°,  structure  XVIII),  202  (5ai03b2^2bi0la22,  structure  XIX), 

201  (5ai23b2°2bi^1a2°,  structure  XX),  201  (5ai°3b222bi^1a2°,  structure  XXI), 
201  (5ai03b2°2bi  11322,  structure  XXII),  2A2  (5ai23b202bi0la2i,  structure 
XXIII),  2A2  (5ai03b222bi0la2i.  structure  XXIV),  2A2  (5ai03b2°2bi2ia2i, 


structure  XXV).  Quartet  coupling  of  the  spins  leads  to  four  states: 
(5ai’3b2^2biHa2°,  structure  XXVI),  (5a1^3b2^2b101a2^  structure  XXVII), 
^02  (5ai^3b2°2biMa2^  structure  XXVIII),  {5a103b2^2b1Ma2^  structure 
XXIX).  The  optimal  geometries  and  relative  energies  of  all  of  these  structures  at 
the  SCF/6-31 1+G*  level  are  presented  in  Fig.  1.  According  to  the  SCF/6- 
31 1  -t-G*  calculations,  XXI  is  the  global  minimum  but  a  few  other  structures  lie 
within  a  few  kcal/mol. 

At  SCF/6-31 1+G*  level,  XV,  XVI,  XVII,  XX,  XXI,  XXII,  XXV,  XXVI,  XXVII 
and  XXVIII  are  local  minima,  XIV,  XIX,  XXIII,  XXIV,  XXIX  are  first  order  saddle 
points,  and  XVIII  is  second  order  saddle  point.  Our  search  for  the  global 
minimum  was  subsequently  expanded  to  0$  symmetry  structures  and  four 
additional  doublet  (2A'  and  ^A")  and  quartet  (^A*  and  M")  structures  were 
studied  at  the  SCF/6-31 1 +G*  level.  The  2a"  state  collapsed  to  XX  or  XXI 
depending  on  the  starting  geometry.  The  2A’  state  has  a  local  minimum  with  the 
structure  XXX  (Fig.  1),  but  this  lies  19.3  kcal/mol  higher  in  the  energy  than 
global  minimum  at  this  level.  The  ^A"  state  converges  to  the  lowest  quartet 
structure  XXVI  and  ^A'  state  has  a  local  minimum  with  the  structure  XXXI  (Fig. 
1)18.2  kcal/mol  higher  in  energy  than  global  minimum. 

Several  of  the  lowest  energy  strucUjres  (XIV,  XVII,  XX,  XXI,  XXIV,  XXVIII 
and  XXIX)  identified  as  detailed  above,  have  been  optimized  at  the  MP2(full)/6' 
31 1 +G*  level  (Fig.  1)  and  the  relative  energies  of  the  rest  of  the  structures  have 
been  calculated  at  MP2(full)/6-311+G*  level  using  UHF/6-311+G*  geometry. 
According  to  our  calculations,  XXI  and  XXIV  has  the  lowest  energies  at  MP2 
level.  The  XXIV  (C2v .  ^Aa)  structure  is  predicted  to  be  the  lowest,  with  the  XXI 
(C2v .  ^Bi)  only  0.05  kcal/mol  (MP2(full)/6-31UG*)  higher  in  energy.  However, 
only  the  XXI  structure  is  a  true  minimum  (its  geometry  is  shown  in  Fig.  1  and  its 
simulated  IR-spectrum  is  presented  in  Fig.  2)  while  the  XXIV  structure  is  a 


saddle  point,  whose  eigenvector  of  imaginary  frequency  produces  an  in-plane 
distortion.  Further  geometry  optimization  of  this  SiaN  structure  within  Cs 
symmetry  was  not  successful  due  to  poor  SCF  convergence.  However  we  do 
not  expect  structures  with  Cs  symmetry  will  differ  much  in  energy  from  the 
energy  of  XXIV. 

Both  neutral  SiaN  and  SiaN  species  have  two  isomers  that  differ  by  only 
a  few  kcal/mol  which  is  different  from  SiaC,  SiaC  and  SiaO,  SiaO,  where  the 
ground  state  structure  is  well  separated  in  energy  from  other  structures. 

Our  computed  adiabatic  ionization  energy  of  SiaN  is  6.3  eV  and  the  most 
energetically  favored  dissociation  energy  is  73  kcal/mol  into  SiaN  Si,  which  is 
much  lower  than  the  dissociation  energy  of  SiaN  into  SiN  •¥  Si  (see  Table  I).  All 
these  number  should  be  considered  as  a  preliminary  estimate. 

Based  on  the  energy  ordering  of  the  neutral  SiaN  molecule  at 
MP2(full)/6-31 1+G*  level  (Fig.  1),  we  studied  two  low-lying  anionic  states  of 
SiaN*:  structure  XXXIII  (Cav,  ^Ai  :1ai2iba22ai23ai2ibi22ba24ai2  3ba2laa2) 
and  structure  XXXIII  (Cav.  ^Ai  :1ai2iba22ai23ai2  Ibi22ba24ai23ba22bi2) 
(geometries  are  presented  in  Fig.  1  ).  However  both  of  these  structures  are 
found  to  not  correspond  to  local  minima.  The  XXXII  structure  has  two  imaginary 
frequencies  and  XXXIII  has  one  imaginary  frequency.  The  latter  structure  is 
lower  in  energy  and  its  eigenvector  of  imaginary  frequency  produces  an  in- 
plane  distortion.  Geometry  optimization  of  the  SiaN*  starting  within  Cs  symmetry 
produces  the  global  minimum  structure  XXXIV  (Cs,  ^A')  structure  whose 
geometry  is  presented  in  Fig.  1  and  the  simulated  IR-spectrum  Is  shown  in  Fig. 
2.  The  final  XXXIV  geometry  of  SiaN*  can  be  considered  as  a  chain  SiiNSi2Si2 
structure.  The  barrier  (XXXIII  structure)  on  the-intramolecular  rearrangement  of 
the  SiaN*  anion  from  one  XXXIV  structure  into  another  is  not  high  (8  kcal/mol  at 
MP2(full)/6-311+G*).  The  structure  XXXV  (C2v.  ^Ai,  1ai22ai2lb223ai2 


1bi24ai22b225ai22bi2)  is  a  saddle  point  with  an  even  lower  barrier  (7.5 
kcal/mol)  for  the  intramolecular  rearrangement  of  the  SiaN*  anion  from  one 
XXXIV  structure  into  another. 

Our  computed  adiabatic  EAa(Si3N)  is  predicted  to  be  1 .4  eV.  The 
dissociation  energies  (at  QCISD(T)  see  Table  1)  for  SiaN*  into  Si-  +  Si2N  is  78 
kcal/mol. 

We  also  optimized  Si2‘*‘,  Si2.  Si2'.  Sia'*'.  Sia  and  Sia'  (simulated 
vibrational  spectra  for  Sia'*'.  Sia  and  Sia*  are  shown  in  Fig.  2),  because  we 
needed  to  know  the  energies  of  these  species  at  a  consistent  level  of  accuracy 
for  the  calculations  of  different  dissociation  channels  of  the  SinN  neutral  and 
ionic  species.  All  of  these  silicon  clusteie  have  been  previously  studied 
experimentally27.28  and  theoretically ,29  and  our  results  for  the  silicon  dimer  and 
trimer  agree  well  with  the  previous  calculations  and  available  experimental 
data. 

V.  MASS  SPECTROMETRIC  RESULTS 

In  order  to  obtain  structural  information,  the  SiaN-*-  cations  were 
subjected  to  a  collisional  activation  (CA)  experiment.  The  resulting 
fragmentation  pattern,  obtained  upon  colliding  the  8  keV  cluster  ions  with 
helium  atoms,  is  shown  in  Figure  3. 

The  base  peak  in  the  CA  spectrum  corresponds  to  the  loss  of  a  neutral 
silicon  atom.  This  Si2N'*‘  fragment  (m/z  70)  contributes  to  approximately  96  %  of 
the  total  intensity  of  the  fragmentation  processes.30  In  addition  to  this  greatly 
favored  dissociation  process,  weak  signals  for  the  losses  of  N,  SiN  and  Si2N 
are  observed  (see  the  20-foid  enlarged  region  of  Rgure  3.).  The  ionic 


fragments  correspond  to  Si3+  (m/z  84).  Si2'^  (m/z  56)  and  Si+  (m/z  28), 
respectively.  Furthermore  the  collisional  activation  process  leads  to  a  rather 
small  fraction  of  SiN+  fragments  (m/z  42). 

For  further  insight  into  the  connectivity  of  the  SiaN^  cluster  ion  and  to 
prove  the  existence  of  the  neutral  SisN  molecule  the  cation  was  subjected  to  a 
neutralization-reionization  experiment,  and  the  resulting  NR-spectrum  is  given 
in  Fig.  4. 

The  NR-spectrum  displays  some  distinct  differences  when  compared  to 
the  results  of  the  CA  experiment.  The  presence  of  a  recovery  signal  for  the 
reionized  SiaN-*-  cations  clearly  demonstrates  that  the  neutral  SiaN  cluster  is  a 
viable  species  in  the  gas  phase.  The  overall  fragmentation  pattern  is  similar 
to  that  observed  in  the  CA  process.  The  SiaN-*-  ion  contributes  strongest  to  the 
total  ion  current.  However,  a  considerably  stronger  fragmentation  of  the 
reionized  SiaN+  into  Sia'*’,  SiN+  and  Si-*-  can  be  found.  In  contrast  to  the  CA 
mass  spectrum  (Fig.  3),  no  signal  at  m/z  84  that  would  correspond  to  the  Sia"*-. 
cation  is  observed  at  the  detection  limit  of  the  machine. 

VI.  DISCUSSION 

According  to  our  calculations,  SiaN,  SiaN*,  and  SiaN  form  two  stable 
isomers  (linear  and  cyclic)  that  differ  in  energy  only  by  few  kcal/mol,  so  both 
may  coexist  in  the  gas  phase.  The  SiaN-*-,  SiaN^,  and  SiaN*  clusters  form  only 
one  very  stable  structure,  so  these  species  should  have  a  single  isomer  in  the 
gas  phase. 

The  theoretical  results  predict  the  most-stable  SiaN-^  cluster  to  possess  a 
planar  rhombic  Cav  structure  (MP2(full)/6-31 1+G*).  In  this  Ui  state,  the  three 
silicon  atoms  do  not  have  identical  Si-N  bond  lengths.  Rather,  two  of  the 


silicon  atoms  are  bonded  equivalently  (Re(Si-N)=1 .720  A),  and  the  third  silicon 
atom  is  coordinated  to  the  nitrogen  atom  by  a  considerably  longer  bond  (1 .894 
A).  The  results  of  the  mass  spectroscopic  experiments  can  be  rationalized  by 
comparing  the  relative  abundance's  of  the  fragment  ions  with  the  theoretically 
derived  reaction  energies  (see  Table  I).  The  calculated  values  predict  the 
fragmentation  process  of  Si3N+  into  Si2N+  and  Si  to  be  highly  favored 
energetically.  This  decomposition  pathway  is  favored  over  the  next  lowest  lying 
reaction  channel  (fragmentation  into  SiaN.  and  Si+)  by  45  kcal/mol  (the  AE 
values  are  all  taken  from  the  results  at  QCISD(T)  calculations).  The  p- 
SisN-*-  SiaN-*-  -t-  Si  can  be  viewed  upon  as  a  simple  silicon  atom  evai.-  ation 
from  the  cluster  molecule  SiaN-^.  All  other  fragmentation  channels  were 
calculated  to  be  significantly  higher  in  energy.  These  results  are  in  excellent 
agreement  with  the  experimental  behavior  of  this  ion  as  the  SiaN-*-  ion  accour  s 
for  96  %  of  the  total  ion  current  in  the  CA  process.  Nitrogen  extrusion  from 
Si3N+  to  produce  the  cyclic  Si3+  (Cav  .^Bi)  requires  60  kcal/mol  more  than  the 
SiaN+  +  Si  formation.  The  CA  spectrum  displays  only  a  very  weak  signal  for  this 
process.  In  the  lower  mass  region  of  the  spectrum,  signals  for  Sia"^,  SiN+  and 
Si-^  are  observed.  The  reaction  energies  for  these  processes  were  calculated  to 
lie  92, 143  and  45  kcal/mol  higher  in  energy  than  the  corresponding  lowest 
channel.  The  "strong''  abundance  of  the  Sia'*’  signal  relative  to  Si+  and  SiN+, 
can  be  rationalized  by  an  additional  contribution  of  secondary  fragmentation 
reactions  of  SI3+  and  SiaN+  that  can  form  Sia"^  upon  collisional  activation. 

The  weak  SiN-*-  signal  agrees  well  with  the  energetic  predictions  as  this 
reaction  channel  requires  1 43  kcal/mol  more  than  the  SiaN"*-  formation.  The 
extrusion  of  a  nitrogen  cation  that  leads  to  the-formation  of  the  cyclic  Sia  species 
(Cav  ^Ai),  requires  251  kcal/mol  more  energy  than  process  1  and  cannot  be 


detected  at  all.  Thus,  the  overall  fragmentation  pattern  is  in  satisfying  agreement 
with  the  expectation  based  on  the  theoretical  values. 

The  SiaN  radical  was  found  to  form  two  different  isomers  that  differ  in 
energy  only  by  a  few  kcal/mol.  The  linear  SiNSi  structure  was  predicted  to  be 
the  global  minimum,  whereas  the  cyclic  structure  lies  only  4.5  kcal/mol  higher. 

The  NR  spectrum  gives  rise  to  a  relatively  small  intensity  for  the  survivor 
ion;  this  observation  points  towards  quite  unfavorable  Franck-Condon  factors 
for  the  charge  transfer  processes.30.3i  The  lowest  energy  requiring 
fragmentation  of  SiaN  is,  in  analogy  to  the  cationic  cluster  molecule,  silicon 
atom  extrusion  accompanied  by  formation  of  the  SiaN  radical.  The  reaction 
energy  for  the  formation  of  SiaN  and  Si  was  calculated  to  be  only  73  kcal/mol. 
The  next  accessible  reaction  channel  (SiaN  Sia  SiN)  requires  53  kcal/mol 
more  than  this  process.  Thus,  again  in  agreement  with  the  experimental  results, 
the  reionized  SiaN-^  species  represents  the  major  dissociation  pathway.  In 
contrast  to  the  fragmentation  in  the  CA  process,  the  NR  spectmm  does  not  give 
rise  to  Sia***  ions.  This  observation  points  to  the  fact  that  the  NR  spectrum  may  be 
viewed  as  a  superposition  of  the  two  mass  spectra  of  Si**-  and  SiaN-*-. 

VII.  CONCLUSIONS 

The  cationic  and  neutral  SiaN  cluster  molecules  have  been 
characterized  by  mass-spectrometric  techniques.  Their  fragmentation  behavior 
in  high  energy  collision  experiments  corresponds  well  with  the  theoretical 
predictions.  High-level  ab  initio  calculations  have  been  used  to  describe  the 
energetic  and  structural  features  of  these  new  molecules.  In  addition,  ail 
experimentally  observed  ionic  and  neutral  fragment  molecules  (SiaN,  SiN,  Sia. 


Si2)  have  been  theoretically  characterized  and  their  vibrational  spectra  are 
reported  in  order  to  assist  their  future  infrared  spectroscopic  characterization. 

According  to  our  calculations,  the  SiaN,  and  SiaN  molecules  form  two 
isomeric  structures  that  differ  by  a  few  kcal/mol  which  may  coexist  in  the  gas 
phase,  while  the  other  clusters  (SiaC,  SiaC,  SiaO  and  SiaC)  possess  one 
preferred  structure. 

Upon  completion  of  this  work,  our  colleague,  Prof.  Michael  D.  Morse, 
informed  us  that  he  and  his  co-worker  0.  J.  Brugh  have  determined^^  that 
ground  electronic  state  of  SiaN  has  linear  structure  with  the  ^Tlg  ground 
electronic  state.  Their  preliminary  equilibrium  bond  length  Re(Si-N)  is  1 .639  A  is 
very  close  to  our  value  Re(Si-N)=1.631  A  at  MP2(full)/6-31 1+G*  and  Re(Si- 
N)=1.643  A  at  QCISD/6-311+G*  levels. 
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Figure  captions 

Fig.  1.  Computationally  optimized  geometries  at  MP2(full)/6-31 1+G*  level 
of  the  neutral,  cationic  and  anionic  SinN  and  Sin  (N=2,3)  species. 
Geometries  at  SCF/6-31 1+G*  are  given  in  parentheses.  Bond  lengths 
are  in  A  and  bond  angles  in  degrees.  Relative  energies  (in  kcal/mol)  of 
SiaN,  SiaN*  at  theSCF/6-311+G*  (in  parentheses)  and  MP2(full)/6- 
31 1 +G*  levels  are  given  under  the  structuries. 

Fig.  2.  Simulated  calculationally  obtained  IR-spectra  of  the  neutral,  cationic  and 
anionic  SinN  and  Sin  (Ns2,3)  species  at  MP2(full)/level,  except  structure 
III  of  SiaN  molecule,  where  calculations  were  made  at  QCISD/6-31 1+G* 
level  (see  text). 

Fig.  3  Collisional  activation  (CA)  mass  spectrum  of  SiaN+  (collision  gas  helium; 
80%  transmission  (T) ). 

Fig.  4.  Neutralization-reionization  (NR)  mass  spectrum  of  SiaN+  (xenon,  80%  T 
//  oxygen  80%  T). 
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